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ABSTRACT 

Using large numbers of simulations of the microwave sky, incorporating the Cosmic 
Microwave Background (CMB) and the Sunyaev-Zel'dovich (SZ) effect due to clusters, 
we investigate the statistics of the power spectrum at microwave frequencies between 
spherical multipoles of 1000 and 10000. From these virtual sky maps, we find that 
the spectrum of the SZ effect has a larger standard deviation by a factor of 3 than 
would be expected from purely Gaussian realizations, and has a distribution that 
is significantly skewed towards higher values, especially when small map sizes are 
used. The standard deviation is also increased by around 10 percent compared to 
the trispectrum calculation due to the clustering of galaxy clusters. We also consider 
the effects of including residual point sources and uncertainties in the gas physics. 
This has implications for the excess power measured in the CMB power spectrum by 
the Cosmic Background Imager and BIMA experiments. Our results indicate that the 
observed excess could be explained using a lower value of erg than previously suggested, 
however the effect is not enough to match eg = 0.825. The uncertainties in the gas 
physics could also play a substantial role. We have made our maps of the SZ effect 
available online. 

Key words: galaxies: clusters: general - cosmic microwave background - large-scale 
structure of Universe 



1 INTRODUCTION 

The discovery of the Cosmic Microwave Background (CMB) 
by Penzias & Wilson ( 1965 ) helped provide the founda- 



tions of modern cosmology. The discovery of the primordial 
anisotropies on large scales by COBE ( Smoot et al.|[l992 |, 
followed by the precision measurements by the Wilkinson 



Microwave Anisotropy Probe ( WMAP; sec for example Hin- 
|shaw et al.||2009 I and other experiments, pinned down the 
values of the principle cosmological parameters. With the 
forthcoming measurements by the Planck satellite, the pri- 
mordial anisotropies at large scales will be measured as ac- 
curately as possible ( Planck consortium|2005 |, with cosmic 
variance preventing any more precise measurements of the 
cosmological parameters from these anisotropies. Interest is 
now turning towards measurements of the primordial polar- 
ization and the total power on smaller scales. 

At these smaller scales, "secondary" anisotropies start 
to dominate over the primordial signal. A significant com- 
ponent of this secondary anisotropy will be caused by the 
Sunyaev-Zel'dovich (SZ) effect. This is the inverse Comp- 
ton scattering of photons from the CMB by free electrons 
( Sunyaev & Zel'dovichll 19701 119721 see iBirkinshawl 19991 for 



a review). It is expected to be dominated by groups and 
clusters of galaxies due to the high temperatures of the in- 
tracluster medium. At microwave frequencies (~ 30 GHz) 



the SZ effect will contaminate the power spectrum due to 
the primary CMB anisotropies at spherical multipoles of 
I ~ 1000 and above, making it difficult to extract cosmo- 
logical information from the primary anisotropies at these 
scales; by I ~ 2000 it is expected that it will dominate the 
power spectrum. 

Several experiments have measured the power spectrum 
in this region of multipoles, and some of these show hints 
of excess power over what would be expected from the pri- 
mary anisotropies alone. The Cosmic Background Interfer- 
ometer (CBI; |Mason et al.|2003||Readhead et al.|2004||Siev 



ers et al. 2009), which observed multipoles up to I — 4000 



at 30 GHz, was the first to measure a possible excess, and 
this has been confirmed by subsequent analysis. This was 
followed by the BIMA interferometer ( Dawson et al.|[2006 l 
observing around I — 5000 at 28.5 GHz, and at higher fre- 
quencies by ACBAR, which has weakly detected a possi- 



ble excess at around I = 2500 at 150 GHz (Reichardt et 



al. 20081. Curiously, observations at similar multipoles by 



the Sunyaev-Zel'dovich Array at 30 GHz (SZA; Sharp et al. 
2009[ ) and QUaD at 100 and 150 GHz (Friedman et al.|2009 1 
do not show any excess power. 

The amount of signal from the SZ effect depends on 
the amplitude of the power spectrum of the density fluctu- 
ations in the Universe. This power spectrum is normalized 
by as, which is the variance of the fluctuations on scales of 
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8h 1 Mpc. To explain the excess measured by the CBI, Siev- 



ers et al. (2009) find that erg must be between 0.9 and 1.0. 



However, observations of the primordial CMB anisotropies 
by WMAP in combination with other instruments and meth- 
ods yield a 8 = 0.812 ± 0.026 ( |Komatsu et al.|2009 |. 

An important issue to consider when estimating the 
range of possible values of erg that agree with the data is 
the statistical properties of the SZ effect, which will be 
non-Gaussian. Here we investigate those properties to see 
whether the excess could be the result of measuring part of 
the sky where the power from the SZ effect is higher than 
average, which would bias the measurement of erg to higher 
than the average value. To do this, we use 3x3 degree 
sky maps with a resolution of 18 arcseconds (600 pixels to 
a side), which allow us to accurately probe the multipolc 
range I = 1000 - 10000. 

To ease comparisons with the observations by CBI and 
BIMA of the high multipole excess, we use a canonical fre- 
quency of 30 GHz. We use a ACDM cosmology throughout 
this paper, with fi m = 0.3, Qa ~ 0.7 and fib — 0.05. The 
Hubble constant Ho = 100 h kms _1 Mpc _1 , where h = 0.7. 
We use three values of erg: 0.75, 0.825 and 0.9. These span 
the best-fitting values from WMAP after 1, 3 and 5 years 



2009) 



of observations (Spergel et al. 2003 2007 Komatsu et al 



For each cosmology, we perform 100 Pinocchio 
( Monaco et~aL|2002| |Monaco, Theuns fc Taffonij|2002| |Ta£ 



foni, Monaco & Theuns 2002 1 simulations to simulate the 



growth of large-scale structure in the Universe; these simu- 
lations are much quicker to run than N-body simulations yet 
provide cluster distributions with an appropriate two-point 
correlation function on large scales. The output of these sim- 
ulations is then used to create 1000 light cone cluster cat- 
alogues. Using a simple model of the gas physics, virtual 
sky maps of the SZ effect are created from the cluster cat- 
alogues and combined with maps of the CMB and residual 
foreground point sources; the latter is based on the extrap- 
olated source number counts from a range of surveys. 

We describe the models and methods used to create the 
virtual sky maps in Section 2. The statistics of these maps 
and their dependency on the cluster model parameters are 
analysed in Section 3. Section 4 discusses the effects of point 
sources, followed by the implications for the high multipole 
excess in Section 5. We conclude in Section 6. 
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Figure 1. The theoretical (red line) and binned (blue points) 
power spectrum of the lensed CMB between I =1000 and 10000, 
averaged over 1000 maps and binned with Al = 1000. The binned 
power spectrum agrees well with the input power spectrum. 
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Figure 2. The statistics of the CMB realizations in the multipole 
bin 2000-3000. The solid red line is the best fit Gaussian distri- 
bution, whilst the dotted blue line is the distribution for this bin 
based on the analytic formula for the variance (equation [2j. The 
distribution is well fit by a Gaussian distribution with a slightly 
higher standard deviation than the analytical fomula. 



2 VIRTUAL SKY MODEL 

The virtual sky maps are built up from a series of compo- 
nents. First, maps of the CMB are generated using a Gaus- 
sian random field. We then use Pinocchio simulations to 
calculate the positions and masses of galaxy clusters, and for 
each cluster we use a simple, spherically symmetric model 
for both the SZ effect and the surface matter density, the 
latter of which is used for distributing point sources (see 
Section B. Full details of each stage are given below. 

The effects of other foregrounds, including the galaxy, 
are not included as the observations focus on areas of the sky 
where galactic foregrounds are at a minimum. Additionally, 
these foregrounds are expected to be fairly smooth on the 
scales under consideration here. 

Similar work creating maps of the SZ effect has been 



carried out by e.g. Geisbiisch & Hobson (2006) 

aT 



(2006) and Sehgal et al. (2007 1; in particular, Holder, 



Schafcr et 



McCarthy & Babul (2007 1 also use Pinocchio simulations 



to create their cluster catalogues. 



2.1 Cosmic Microwave Background anisotropies 

Maps of the small-scale CMB anisotropies can be computed 
using the standard assumption that the CMB is a Gaussian 
random field and hence is completely described by its an- 
gular power spectrum, Ci. We generate the expected power 
spectrum of the gravitationally lensed CMB from CAME^] 
(Lewis, Challinor & Lasenby 20001. It is known that the 
gravitational lensing component is non-Gaussian (see for 



1 November 2008 version; available from http://www.camb.info 
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example Zaldarriaga 2000), however as this component is 
sub-dominant we assume that it is Gaussianly distributed. 

When calculating the binned power spectrum, we as- 
sume that 1(1 + l)C;/(27r) is flat over a region of size Al 
such that the power within a multipole bin centred on U 
can be calculated by 



B r = 



1 

a7 



E 

;ebin 



1(1 + l)Ci 
2tt : 



(1) 



where the summation is over Ti — |AZ ^ I < h + ^Al. Values 
of Ci that have not been sampled due to the finite grid are 
assumed to have the value of the closest Ci that has been 
sampled. This is obviously imperfect where the Ci's have a 
large gradient or are ill- or irregularly-sampled, but is likely 
to be sufficient for the purposes at hand. 

The expected mean in each multipole bin can be cal- 
culated from the input power spectrum, and the expected 
variance is calculated by 



SBi 



1 

~Al 



\ 



Zebin 



l(l + l)Ci 
2tt 



(2/ + l)/sky • 



(2) 



where / s ky is the fraction of the sky, Al is the width of the 
multipole bin (here, Al = 1000) and the sum is over all mul- 



tipoles within the bin. The quantity SBi 



i: 



1/2 .-1/2 



sky 



sometimes called the cosmic variance since it represents the 
limit on how well one can possibly measure Bi as imposed 
by the Gaussian statistics. The aim of this paper is to under- 
stand this issue in the case of the non-Gaussian anisotropies 
created by the SZ effect. 

Fig. [l] shows the CMB power spectrum using bins of 
1000 multipoles, and Fig. [2] shows the statistics within the 
multipole bin 2000-3000. The statistics are consistent with a 
Gaussian distribution with a slightly larger standard devia- 
tion than expected from cosmic variance due to inefficiencies 
in the sampling of the spectrum at different multipoles. The 
distribution is also slightly skewed, as the finite sampling 
leads to a \ 2 distribution rather than a Gaussian one. 



2.2 Cluster catalogues 

In order to create a map of the SZ effect, the positions and 
properties of the clusters need to be determined. We use 
the PINpointing Orbit-Crossin g Collapsed Hierarchical Ob- 
jects (Pinocchio) algorithm^ ( Monaco et al.|2002 Monaco, 



Theuns fc Taffoni|2002[ |Taffoni, Monaco fc Theuns|2002[ ) to 
generate mock cluster catalogues. 

Pinocchio uses Lagrangian perturbation theory to pre- 
dict the collapse of matter into haloes. Galaxy cluster po- 
sitions and masses generated by Pinocchio will be associ- 
ated with the large-scale cosmological structure, and hence 
will have a non-zero two-point correlation function (they 
will be clustered), whereas distributions created from ana- 
lytical mass functions (for example Sheth fc Tormen||1999 



henceforth ST99) with random cluster positions will be Pois- 
son distributed with zero two-point correlations. Pinocchio 
does not simulate the structure within haloes, thus offering 



2 Version 2.1.2 beta, available from 

http : //adlibitum. oats . inaf . it/monaco/Homepage/Pinocchio/ 
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Figure 3. The ratio of the PINOCCHIO mass function to ST99 
for different resolutions of the same PINOCCHIO simulation. The 
resolutions are 2.5 Mpc (blue short dashed line), 1.66 Mpc (black 
dotted line) and 1.25 Mpc (red solid line; our adopted resolu- 
tion). The latter two are consistent with each other down to 
A/ vir ~ lO 13 A/ . At Mvir = 10 13 Af Q the mass function from 
the highest resolution simulation is in good agreement with that 
from ST99 for a resolution of 1.25 Mpc. At the highest masses 
the scatter is caused by small number statistics, but a systematic 
trend downwards is still present. 



a significantly faster method of creating cluster catalogues 
compared with N-body techniques. In contrast to N-body 
simulations, requiring many thousands of CPU-hours to run, 
a reasonable resolution Pinocchio run takes less than a day 
on a single CPU. 

We have run 100 Pinocchio simulations for each value 
of erg, each simulating comoving cubes with sides of 500 Mpc 
with a cell size of 1.25 Mpc (thus 400 cells to a side). Halo 
catalogues were created at 50 redshifts between z = 5 and 
z = 0, the centres of which are spaced 160 Mpc apart in 
comoving distance. The initial power spectrum used to gen- 
erate the distribution of matter in the simulation is from 



Efstathiou, Bond & White (1992) 



A set of 1000 lightcone catalogues for each value of as 
were generated from the simulation output. For each of the 
50 redshifts, a halo catalogue is randomly selected from the 
100 simulations. This halo catalogue is then translated by 
random amounts, independently in all three dimensions us- 
ing a periodic box, as well as being rotated by a randomly 
chosen integer multiple of 90 degrees independently along 
each of the three axes. A 160 Mpc deep slice of the sim- 
ulation is then taken, and clusters within a 3 by 3 degree 
block centred on the middle of the slice are included in the 
lightcone. The limit of 3 degrees to a side is imposed by the 
size of the simulation at the highest redshift; the box size 
subtends an angle of 3 degrees at z = 5. 

The choice of simulation resolution means that the sim- 
ulations have a high degree of completeness down to M V ir = 
1O 13 M0; that is, the comoving number density of clusters 
n(M) = dn/dM agrees between the Pinocchio simulations 
and that expected from ST99, as shown in Fig. [3] However, 
Pinocchio systematically underestimates the number of the 
largest mass clusters when compared to the analytical pre- 
diction of ST99; this becomes more pronounced at higher 
redshifts (see Fig. [41. This is unexpected since Pinocchio 
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Figure 4. The comoving differential number count of clusters 
greater than a given mass analytically from ST99 (dashed lines), 
from the mean mass function of all PlNOCCHlO catalogues with 
<rg = 0.9 (red +) and from the corrected PlNOCCHlO catalogues 
(blue crosses). The number counts for four redshifts are given; 
from the top down these are z=0, 1, 2 and 3. The correction 
brings the two mass functions into agreement 



is designed to reproduce ST99, and causes a significant dis- 
crepancy in the amplitude of the power spectrum of the SZ 
effect compared with maps created from random cluster cat- 
alogues generated from the ST99 number count predictions. 
We assume that the ST99 mass function is the more ac- 



curate of the two, although Warren et al. (2006) show that 



this may itself be systematically low at the highest masses 
compared to N-body simulations. We calculate a mass cor- 
rection table for each cosmology such that the number of 
clusters greater than a given mass at each redshift (calcu- 
lated from the average of all of the simulations with the 
same cosmology) matches the predicted value from ST99. 
The appropriate correction is then applied to all clusters in 
the lightcone prior to the creation of the maps. 

The effect of this correction is shown in Fig. [4] which 
gives the differential number count of clusters greater than 
a given mass before and after the correction has been ap- 
plied. The correction brings the PlNOCCHlO spectrum into 
agreement with that generated from a set of 1000 cluster cat- 
alogues with mass and redshift values drawn from the ST99 
mass function and randomly determined positions. The av- 
erage mass and redshift distributions of the corrected cata- 
logues for all three cosmologies are given in Fig. [5] 



2.3 Cluster model 

Although clusters are complicated in the optical regime, for 
the purpose of modeling the SZ effect they can be repre- 
sented very simply by assuming that the clusters are relaxed, 
virialised objects. Although clusters undergo major mergers, 
these are relatively infrequent, occurring in around 10 to 20 



percent of the population depending on the redshift ( Kay et 
a^[2007t 



The cluster model we utilize is spherically symmetric 
and consists of two components: the dark matter and the 
gas. The former is used to calculate the surface mass density 
of the cluster, used for the inclusion of point sources, and the 
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Figure 5. The distribution of all clusters as a function of redshift 
(top) and mass (bottom) within the modified lightcone catalogues 
above 10 13 Mq, for each of the three values of erg (from the top, 
0.9 with black dots, 0.825 with the red solid line and 0.7 with blue 
large dots). Values from ST99 would overlay the points plotted; 
they are essentially the same (by construction). 



latter gives the thermal SZ effect. A kinetic SZ effect also 
exists, which is created by the relative motion of the clusters 
to the CMB, but this is likely to be negligible compared to 
the thermal SZ effect at the frequencies discussed in this 
paper and hence is not considered here. 

The virial radius, r v i r , is used to define the size of a 
cluster. This can be calculated by using the spherical infall 
model ( |Gunn fc Gott|[l972[ ), which gives 



f vir 



3M v ir 



4TYA c (z)p crit (z) 



1/3 



(3) 



where M v i r is the virial mass of the cluster, p C rit(z) = 
3H 2 (z)/ (8nG) is the critical density for the Universe to be 
flat and H(z) is the Hubble parameter, given by the Fried- 
mann equation 



E\z) 



H 2 (z) 



On(l + 



(4) 



for a flat universe (which we will assume throughout). The 
quantity A c is the mean matter density within the virial 
radius of the cluster in units of the critical density, 
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A c (z) = 



pcluster (^) 
Pcrit(z) 



(5) 



The value for A c depends on the cosmological parameters; 
in an Einstein-de Sitter universe it is exactly 187T 2 » 178. 
For our adopted cosmology, we use the fit given by |Bryan fc| 
Narmanl (11998b, 



18tt 2 + 82(Qm(z) 



(6) 



where Q m (z) = Q, m (l + z ) 3 /E 2 ( z ) is the matter density at 
redshift z. 

Finally, the cluster radius needs to be converted to an 
angular size so that the cluster can be projected on to a 
virtual sky map. This is done by dividing the virial radius 
by the angular diameter distance, (9 v i r = r^jd^z) in which 



(1a(z) = 



:flK-\z')dz' 
(1+*) ' 



(7) 



2.3.1 Dark matter 



We calculate the surface mass density of each galaxy cluster 
for the purpose of distributing point sources on to the map 
(see Section [2J. For the dark matter density profile, we use 
the NFW profile, 



Pdm 



Pcrit Ac 
i 1 fl+ & ' 



(8) 



This was found by fitting profiles to dark matter haloes from 



N-body/hydrodynamical simulations run by Navarro, Frenk 
|fc White| ( [T995| |l996"l [1997] ), and is now the standard pro- 
file for modeling dark matter haloes. The scale radius for 
the NFW profile is r s — r vlr /cDM, where cdm is the concen- 
tration parameter. Values for the concentration parameter 
depend weakly on the cluster mass with a certain amount of 
scatter (see for example |Neto et al.|2007||Duffy et al.|2008 |. 
However, as we are only using the surface mass density maps 
indirectly these effects are unimportant here, so we assume 
a constant value, cdm = 5, which is appropriate for clusters. 

The surface mass density at distance ip s = 9/9 B from 
the centre of the cluster, where 9 S = r s /dA and 9 is the 
angular distance from any given position in the sky to the 
centre of the cluster, is given by 



(9) 



In this, the central surface mass density of a cluster can be 
calculated using (see for example Lokas fc Mamon|2001 1 

M vir 



So 



27rr2(]n(l + c)-(c/(l + c)))' 



(10) 



and the NFW profile can be projected on to a 2D sky to 
give ( Bartelmann|199"6 1 



Cdm 



• arctan 



<ps-i 



• arctan 



Ifs > 1, 

<P.<1, 

lfi s = 1. 



2.3.2 SZ effect 

Our gas model consists of two parts: a profile, and a nor- 
malization for the level of the SZ effect, which is dependent 



on the mass and redshift of the cluster. The method is anal- 
ogous to that used for creating maps of the dark matter 
distribution. 



We use the isothermal /3-model (Cavaliere & Fusco- 
Femiano|1976 I for the cluster profile, 



= i + 



(12) 



in which /3 is a dimensionless parameter that measures the 



outer slope of the profile, r c 



/c gas is the core radius 



that describes the turn-over point between the core and the 
power-law slope. The concentration of the SZ effect around 
the cluster centre is controlled by c gas . Fiducial values of 
Cgas = 10 and f3 — 2/3 are used, following Battye & Weller 
(2003). We truncate the profile at the virial radius, r v i r , 
which prevents potential divergence in the gas mass (see for 
example Birkinshaw 1999 1. This may however underestimate 



the total SZ effect from an individual cluster, where infalling 
gas is shocked outside of the virial radius and prevented 
from falling inwards (see for example Evrard 1990; Kocsis, 
Haiman fc Frei|2005[ ) 



The change in the temperature of the CMB from the 
SZ effect can be calculated using 



AT(v? c ) = y(ip c )g(x)T C MB 



(13) 



where g(x) — (xj tanh(s/2)) — 4, the dimensionless fre- 
quency x — hpu/kBTcMB, is Boltzmann's constant, hp 
is Planck's constant and v is the frequency of interest. The 
present-day temperature of the CMB as measured by the 
FIRAS i nstrument on COBE is T C mb = 2.728K ( |Fixsen et 
al.|1996 1. The SZ effect along the line of sight from a single 



cluster depends on the integrated SZ effect from the cluster, 
Y(M,z), and the cluster profile C(ip c ) as 



V(<Pc) = 



2tt J^' a<P'cWcd<p' c 



(14) 



where <p c = 8/9 c is the distance from the centre of the cluster 
in units of 9 C = r c /d\. 

The integrated Y parameter is a measure of the total 
power of the SZ effect and is used as the normalization of 
our cluster profile. It is directly related to the total thermal 
energy of the gas. We assume a power-law relation between 
Y and the cluster mass, 



Y 



Y.h- 1 ( M vir \ ( A c (z) 
d\{z) I lO^-iMe J I A c (0) 



1/3 



E\z)\ (1 + *)°,(15) 



where Y* , 7 and a fix the normalisation, slope and red- 
shift evolution respectively. For our fiducial values we adopt 
Y* = 2 x 10" 6 Mpc 2 ; the isothermal value of 7 = 5/3 and 
a = such that the clusters are approximately self-similar. 
These choices are motivated by the results from recent cos- 
mological simulations, which demonstrate that the relation- 
ship between Y and M is a power-law with small intrin- 
sic scatter, being close to that predicted by the self-similar 
model at all redshifts, and relatively insensitive to the effects 
of cooling and heating of the intracluster medium (see, for 
example, |da Silva et al.|2004| |Motl et al.|2005} |Nagai|2006l ) . 

Using the recent Millennium Gas simulations, Kay et al. 
(in prep.) obtain best-fitting values of Y* = 2.3 x 10~ 6 Mpc 2 
and 7 = 1.64 for a non-radiative simulation at z — 0, and 
Y» = 1.9 x 10 -6 Mpc 2 and 7 = 1.76 for a simulation where 
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Figure 6. A realization from the as = 0.825 cosmology using a clustered PlNOCCHlO lightcone with mass corrections. The maps are lxl 
degrees, with a resolution of 6 arcseconds. From left to right the maps are of the temperature decrement from the SZ effect, the CMB 
and the two combined. White represents colder areas, with black representing hotter areas. The colours range from to -200 fiK for the 
first image and 200 to -200 fiK for the other two. 



the gas was preheated and allowed to cool radiatively. Both 
fits were applied to clusters with M v i r > W 14 Ii~ 1 Mq, the 
objects that dominate the power spectrum over the range of 
multipoles of interest in this paper. These values agree well 
with the fiducial values above. 

Finally, the j3 profile can be projected from three di- 



mensions to a 2D sky to become ( Battye & Wcller 2003 ) 



1 3/3 

Csz(^ c ) = + 2 



1/2 



J[c,(3] 



where 

J[a,b]= f ( 
Jo 



1 + x 2 ) 2 dx. 



(16) 



(17) 



2.4 Realization properties 

Fig. [6] shows lxl degree maps of the SZ effect and the 
CMB separately and combined for a single realization from 
the 08 = 0.825 cosmology. In the SZ map, an M v i r = 
5.4 x 1O 14 M0 galaxy cluster at a redshift of z = 0.13 lies 
in the top right, with a virial radius of 14 arcminutes. The 
SZ map contains 1,933 galaxy groups and clusters in total, 
with a maximum decrement in a single pixel of 450/xK and 
an average of 8.5pK per pixel across the map. 

For the 3x3 realizations with og = 0.75 and a resolu- 
tion of 18 arcseconds, there are 13470 ± 300 galaxy groups 
and clusters per realization, where the standard deviation 
is from the scatter between the realizations. These objects 
cause an average decrement per pixel of 4.3±0.5/iK in each 
map (again with a standard deviation from the scatter). For 
o"8 = 0.825 this increases to 18040 ± 350 galaxy groups and 
clusters and an average of 6.5 ± 1.3/iK, and for as ~ 0.9 this 
further increases to 23000T400 galaxy groups and clusters 
and an average of 9.3 ± 1.3/iK. The increase in the number 
of objects, and hence the total amount of power from the 
SZ effect, is due to the increase in the clustering of matter 
as as increases. 

To confirm that the power spectra calculated from these 
realizations are converged such that no significant amount 
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Figure 7. Top panel: The mean power spectra from the SZ effect 
in the erg = 0.825 cosmology imposing various minimum mass 
limits: from the bottom up, 1 X 10 15 , 5 X 10 14 , 2 X 10 14 , 1 X 10 14 , 
5 X 10 13 , 2 X 10 13 and 1 X 1O 13 M . Bottom panel: The same, 
but imposing a maximum redshift cutoff of (from the bottom up) 
z = 0.5, 1.0, 1.5, 2.0, 2.5, 3.0, 3.5, 4.0 and 4.5. These results show 
that the spectrum is approximately converged for M > 10 13 Afo 
and z < 3.5. 
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Figure 8. The power spectra expected from the SZ effect as 
calculated using the analytical formula and the modified PlNOC- 
CHIO realizations for each cosmology. The blue dashed line shows 
<rg = 0.75, the red solid line 0.825 and the black dotted line 0.9. 
The higher values of the pairs are from the theoretical spectrum, 
the lower from the realizations; there is consistency between the 
shapes of the curves and the analytical formula is 2 — 4 per cent 
higher than the realizations. 



of power comes from clusters with lower masses or higher 
redshifts, the mean spectra from realizations with various 
minimum mass and maximum redshift cutoffs imposed are 
shown in Fig.[7]for a s = 0.825. The spectra are converged to 
within a few /iK 2 over the multipole range of interest when 
clusters down to 10 13 Mq and out to a redshift z — 3.5 are 
included (the effects of reducing the maximum mass cutoff 
will be discussed in Section [H]). We use 10 13 M Q and z = 4.5 
as the fiducial limits for the realizations shown in the rest of 
this paper. Although the average SZ effect across the map 
is converged by the same redshift, it is not converged at a 
minimum mass of 1O 13 M0; lower mass objects than consid- 
ered here will significantly contribute towards the average 
SZ effect. 

As an additional check, we compare the ratio of the 
mean power spectra between realizations with different val- 
ues of as- Komatsu & Seljak (2002 1 found that Ci scales 



as <Tg , where a ss 7. Fitting for a for our realizations and 
averaging over the multipole bins, we find that between the 
as = 0.75 and 0.825 realizations, a ~ 7.1. Between as = 0.90 
and 0.825 the value is a w 6.9. Thus there is fairly good 
agreement, although the relationship appears not to be a 
perfect power law over this range of as- 

The mean power spectrum can also be evaluated di- 
rectly from the number density and cluster profile using the 
halo formalism (see for example |Komatsu &: Seljak||2002| ) . 
The details of this calculation are outlined in the Appendix, 
and Fig. [8] compares the values computed using this method 
with those from the realizations for the three different cos- 
mologies. The analytical method predicts ~ 2 — 4 per cent 
more power than the realizations, which is likely due to in- 
efficiencies in the creation of the maps, but it appears that 
there is broad consistency here. 



We have made our maps of the SZ effect with as = 0.825 
available online in FITS format^] 



3 POWER SPECTRUM STATISTICS 

Using the virtual sky maps described in Section 2, the binned 
power spectra from the SZ effect are calculated using bins 
of AZ = 1000. We analyse the histograms of these power 
spectra and investigate the normalized skew as well as the 
mean and standard deviation. Higher orders than the skew- 
ness, such as the kurtosis, tend to be significantly affected 
by outlying data points and hence are not considered here. 

Similar work investigating the statistics of SZ effect 
realizations has been carried out by |Zhang, Pen fc Wang| 
( |2002[ > and |Shaw etaT] ( |2009 K 



3.1 Fiducial results 

The resulting histograms for the combined 3x3 degree CMB 
and SZ effect maps with as = 0.825 between I — 1000 and 
9000 are shown in Fig. [9] To quantify the distributions, we 
compute the mean Bi for each bin i, as well as the stan- 
dard deviation SBi and normalized skewness Si. The latter 
is defined by 



N ^ 



B? - Bi 



(18) 



where B™ is the value of the binned 1(1 + l)C'i/(2n) from 
the nth realization, Bi is the mean of those values, N is the 
number of realizations and the sum is over all realizations. 
The computed values for all three values of as, for both the 
power spectrum of the SZ effect on its own and that from 
the combined CMB and SZ effect, are given in Table [1] 

The SZ effect significantly increases the standard devi- 
ation of the realization statistics for all but the lowest mul- 
tipoles. This is shown graphically in Fig. |10| where the ratio 
of the standard deviation to the mean is given for the CMB 
and SZ effects separately, and then combined. The SZ effect 
has a much higher ratio than for the CMB on its own. When 
combined with the CMB, the realizations have a low stan- 
dard deviation at the lowest multipoles, where the CMB is 
dominant, but quickly return to higher values at the higher 
multipoles where the CMB has decreased in power due to 
Silk damping and the SZ effect is closer to its peak power. 
The point at which the transition occurs is slightly different 
for the three values of as; this is due to the different power 
levels from the SZ effect, with the as = 0.9 cosmology having 
more power and hence transitioning at lower multipoles. 

Fig. [TT] shows the equivalent histograms to Fig. [9] but 
at 150 GHz. The SZ effect at this frequency is weaker by a 
factor of 2 than at 30 GHz, reducing the power spectrum 
by a factor of 4. As a result, the CMB dominates the power 
spectrum to higher multipoles and the increase in standard 
deviation due to the SZ effect does not become important 
until I ~ 4000. The highest multipoles considered here are 
not greatly affected by this change in frequency as these 
remain dominated by the SZ effect. 
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Figure 9. The statistics of the combined CMB and SZ effect power spectra for the <rg = 0.825 cosmology for a thousand 3x3 degree 
maps. In each histogram, the x-axis is the binned power Bi in the maps in /xK 2 , and the y-axis is the probability density of that power in 
a map. The histograms on the top row show bins of 1000 multipoles between I =1000 and 4000; the second row 5000-9000. A Gaussian 
curve generated from the mean and standard deviation of the distribution is shown by the dotted blue line. The red solid line shows the 
equivalent for a set of 1000 realizations generated using the mean power spectrum from the combination of the CMB and SZ effect but 
assuming that all the power is Gaussianly distributed. This agrees at the lowest multipoles, but quickly becomes too narrow, indicative 
of the increased variance due to the SZ effect. 
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Table 1. The statistics of the SZ effect (top section) and combined CMB and SZ effect (bottom section), generated using the 3x3 degree 
clustered galaxy cluster catalogues from the modified PlNOCCHlO simulations. The mean (Bi) and the standard deviation (SBi) within 
bin i are in units of /iK 2 ; the skew (s) is dimensionless. The ratio of the standard deviation to that from the Gaussianly distributed 
simulations with the same power spectrum is given by SBi/SBi^Q. There is a considerable increase in the standard deviation from the 
SZ effect compared to Gaussian distributions, and the distribution is also skewed. 



To quantify more precisely this increase in standard de- 
viation, we create realizations with the same mean power 
spectrum but Gaussianly distributing all of the power in the 
map, as if all of the power is due to the CMB. The Gaus- 
sian fit to these distributions is shown by the red solid lines 
in Fig. [9] and |11| At the lowest multipoles, where the CMB 



is dominant, the statistics from this method are in good 
agreement with those from the realizations. At higher mul- 
tipoles, as expected, the standard deviation is much lower 
for Gaussian statistics than for the SZ effect. The ratios of 
the standard deviation from the two methods are given in 
Table [T] the standard deviation is much smaller, typically 
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Figure 11. As Fig. [9] but for 150 GHz. The CMB dominates to higher multipolcs than at 30 GHz., such that the broadening of the 
distributions due to the SZ effect is less important at intermediate multipoles. 
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Figure 10. The ratio of the standard deviation to the mean for 
the CMB only (solid black line), the SZ effect only (top three 
lines) and the combined CMB and SZ effect (central lines) for 
the three different values of erg (0.75: blue dotted line, 0.825: red 
solid line, 0.9: black double-dotted line) from the 3x3 degree 
maps. The values for the CMB scale as Z -1 / 2 , as per equation]^ 
whereas those from the SZ effect scale as l~ 1 but with a higher 
amplitude. 



Figure 12. The skewness s for the 3x3 degree maps for the three 
values of erg: 0.75 is shown by the blue dotted line, 0.825 by the 
red solid line and 0.9 by the black double-dotted line. The top set 
of lines are for the SZ only; the bottom set are when combined 
with the CMB. There is little difference between the two above 
I ~ 4000. At the lowest multipoles, the SZ effect significantly 
skews the distribution, and this skewness remains important at 
intermediate to high multipoles. 



by a factor of 3 for I > 3000, when compared with the non- 
Gaussian realizations. 

In Section |2.4| the mean power spectra were found to 
scale approximately as Ci oc <7g , where a ~ 7, as per |Ko-| 
|matsu fc Seljak"| ( |2002[ ) . Assuming that the standard de- 
viation varies with ag in the same fashion, we find that 
q sd ss 6.0 between erg = 0.75 and 0.825, and ~ 6.8 be- 
tween 0.9 and 0.825. Thus, there is a strong dependence on 
<rg for the standard deviation, although not quite as strong 
as for the mean power spectrum. The curve for erg = 0.75 
is slightly higher than for the other two values of ag in Fig. 
|10| this is an illustration of the non- universality of qsd- 



In Section |2.4[ the mean theoretical spectrum calcu- 
lated by an analytical formula was compared to the mean 
from the realizations and was found to agree well. The stan- 
dard deviation can also be computed; see the Appendix for 
the method. The standard deviation calculated solely from 
the Gaussian components differs from that calculated from 
the realizations by a factor between 3 and 7 depending on 
the multipole; when the angular trispectrum component is 
added the two agree to within 20 per cent when compared 
with the results from the PlNOCCHlO realizations presented 
in Fig. |10| The origin of the remaining discrepancy will be 
discussed in section [3~4l 
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Figure 13. The correlation matrix for 1000 realizations of the CMB (left), SZ effect (middle) and the combined CMB and SZ effect 
(right), using 3x3 degree maps with bins of Al = 1000. The bin numbers are given in the x and y axis (for example, is the multipole 
range — 1000). The CMB has negligible off-diagonal terms, however these terms are no longer negligible when the SZ effect is accounted 
for. 



We find that the SZ effect results in a positive skew 
in the statistics, resulting from an overabundance of high- 
power realizations when compared with the expectations 
from Gaussian statistics due to the presence of rare, mas- 
sive clusters and associated smaller clusters. This is shown 
in Fig. |12| the lowest multipoles for the SZ effect on its own 
are especially skewed. With the addition of the CMB, the 
low multipole bins become Gaussian as the CMB dominates, 
but some skewness remains at the higher multipoles where 
the SZ effect is dominant. There is a certain amount of noise 
in the measurements of the skew here, however; an increased 
number of realizations would measure this more accurately. 



Zhang fc Sheth (20071 have also found that the SZ effect 



yields a positively skewed distribution using an analytical 
approach for calculating the probability distribution for a 
single multipole. 

3.2 Correlation matrix 

For purely Gaussian maps, for example the CMB, the power 
on different scales is uncorrelated such that the off-diagonal 
terms of the covariance matrix are zero. However, this may 
not be the case for the SZ effect. The covariance matrix can 
be calculated by 



(19) 



This is then normalized to give the correlation matrix by 

(20) 



The correlation matrix is depicted graphically for the 
CMB and SZ effect individually and combined in Fig. |13| 
The matrix is diagonal for the CMB, however the off- 
diagonal terms are significant for the SZ effect, in agreement 



with the predictions of Cooray (20011. This is due to the 



power from individual clusters spanning many thousands of 
multipoles, correlating the multipole bins considered here. 
This effect is reduced for widely separated multipole bins as 
individual clusters are not dominant at all multipoles; for 
example, as shown by Fig. [7] the largest clusters only pro- 
vide the bulk of the SZ effect at the lowest multipoles, and 
smaller (mostly uncorrelated) clusters provide the power at 
the highest multipoles. 



3.3 Dependency on map size 

The statistics of the SZ effect power spectrum will depend 
on the size of the map being considered. We investigate this 
dependency by considering maps of 1 x 1 and 2x2 degrees, in 
addition to the 3x3 degree maps that have been investigated 
thus far. These smaller maps are created by selecting and 
mapping only the clusters from the catalogues that lie within 
these areas. For the lxl degree maps, this provides 9000 
independent realizations, whereas for the 2x2 degree maps 
we continue using 1000 realizations. The mean power spectra 
for these three different map sizes remain the same. 

The histograms for the lxl degree realizations contain- 
ing both the CMB and the SZ effect are shown in Fig. [14] As 
before, for the lowest multipoles the CMB dominates, so the 
histogram for the lowest bin is approximately Gaussian. For 
the higher multipoles, however, the histograms are signifi- 
cantly skewed towards positive values, much more so than 
for the 3x3 degree maps. This qualitatively agrees with 
|Zhang fc Sh eth (2007|, who found that larger map sizes de- 
crease the skewness of the probability distribution due to 
averaging over a larger number of clusters. 

For the CMB, the standard deviation from the realiza- 
tions scales linearly with V/g/y ■ The ratio of the standard 
deviation to the mean for the three maps sizes is shown in 



Fig. 



15 



this scales as 1//. 



1/2 
sky 



for both the SZ effect on its 
own and the combined CMB and SZ effect. This is as pre- 
dicted by the analytical formula. The transition between the 
CMB- and SZ-dominated regimes is clear for all three map 
sizes. The skew for the SZ effect also scales close to l// s ky • 
If the dependence is parameterized as 1/ f^ 2 



'sky 



and a is cal- 
culated independently for each multipole bin, then a w 1.2 
at the lowest multipoles, changing to m 0.9 at the highest 
multipoles. This scaling indicates that the skewness can be 
reduced to < 0.1 for maps larger than 400 square degrees. 

As we are only considering the clusters within the map 
size, we will not be including any effects from clusters that 
are positioned outside the map but extend into the map. 
This was true for the larger maps considered earlier, however 
this may become more important here as large clusters can 
extend up to a quarter of a degree in the maps. As the 
means from the three map sizes differ by less than a percent, 
however, and the standard deviations agree with predictions, 
this appears to be negligible. 
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Figure 14. As Fig. [5] but for 9000 1 x 1 dc gree maps. The first few multipole bins are dominated by the CMB, hence have a distribution 
that is close to Gaussian; the later bins become increasingly more skewed to higher values. Extremal values are more likely than in the 
3x3 degree maps. 
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Figure 15. Left: 8Bi/Bi from the SZ effect on its own (top lines) and for the CMB and SZ effect combined (bottom lines) with erg = 0.825 
for three different map sizes: 3x3 degree (red solid line), 2x2 degree (blue dotted line) and lxl degree (black double-dotted line). The 
values scale as l//g^ 2 ■ Right: As the left plot, but plotting the skewness for the three map sizes; this also scales roughly as l/Z^y ■ 



3.4 Effects of clustering 

To check whether the clustering of the galaxy clusters can 
be responsible for the non-Gaussianity in the statistics, we 
randomize the positions of the clusters on the sky. This ef- 
fectively modifies the angular correlation function to that 
expected of a Poissonian distribution, although this still in- 
cludes the effects of clustering due to fluctuations that are 
larger than the map size. For the 3x3 degree maps, there is 
no significant change in any of the values of the mean, vari- 
ance or normalized skew between the statistics of the ran- 
domized and the clustered maps. However, when the values 
of SBi/Bi are compared to those from the analytical calcu- 
lations (as described in the Appendix), a 10 — 15 per cent 
excess is present at high multipoles 

Simply randomising the positions on the sky effectively 



destroys the clustering between galaxy clusters on scales 
smaller than the fiducial map size of 3 x 3 degrees. However, 
this will not remove the effects of clustering on the mean 
number density of clusters within the maps. To illustrate 
the effect of this, we create sets of 1 x 1 degree maps where 
we randomise the clusters by two methods: first, selecting 
the clusters within a 1 x 1 degree area then randomising 
their positions (analogous to the randomization performed 
on the 3x3 degree maps), and second, randomising the 
whole map and then selecting lxl degree patches. 

Selecting the clusters that lie within the lxl degree 
area and then randomizing their positions provides largely 
the same statistical properties as the clustered realizations, 
as was the case for the 3x3 degree realizations. However, 
if the clusters are randomized prior to selection, then the 
ratio of SBi / Bi is close to agreement with the analytical es- 
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Figure 16. Top: The ratio of 8Bi/Bi from clustered and ran- 
domised lxl degree maps for the three different values of erg 
(the red solid line is erg = 0.825, the blue dotted line 0.9 and 
the black double-dotted line 0.75). Clustering increases the stan- 
dard deviation by ~ 10 per cent. Bottom: The same as the top 
figure, but comparing the clustered realizations to the analytical 
expectation. The increase of ~ 10 per cent remains. 

timate, differing by a few percent. The effect of randomising 
and then selecting the clusters is shown in Fig. |16| which 
gives the ratio of SBi between the clustered and random 
realizations (Bi is the same for the two methods) and be- 
tween the clustered realizations and the analytical values. 
Clustering increases SBi/Bi by ~ 10 per cent at the high- 
est multipoles, where the standard deviation is increased 
as the galaxy clusters group together in certain realizations 
and are absent in others, broadening the standard deviation. 
The lowest multipoles are not well sampled and are highly 
skewed, such that their standard deviation alone does not 
adequately represent the distribution. 

3.5 Upper mass limit 

The power spectrum of the SZ effect measured in different 
parts of the sky is very sensitive to the mass of the largest 
clusters within that field. Fields used for CMB observations 
might be preselected to avoid such objects, hence biasing 
our measurement of the power spectrum. In order to test the 
dependence on this, we apply a series of maximum mass cuts 
to the PlNOCCHlO catalogues; the mean spectra after these 
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Figure 17. Top: The mean power spectra from the SZ effect in 
the <rg = 0.825 cosmology with all clusters (red solid line) and 
after removing all clusters greater than 8x (green large-dashed), 
6x (blue small-dashed), 4x (pink dotted) and 2 X lO 14 h _1 M0 
(turquoise dot-dashed), from the top down. Middle: 8Bi/Bi from 
the SZ effect as a function of the maximum mass in the re- 
alizations in two bins: I = 2000 — 3000 (red solid line) and 
I = 3000 — 4000 (blue dotted line). Bottom: the skewness s of 
the distributions for the same mass cuts as the mean power spec- 
tra. All three of these quantities are significantly affected by the 
largest mass clusters in the realizations. 
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cuts are shown in Fig. |17| Note that the effect of minimum 
mass cuts was discussed in Section \2. 41 

The mean power spectrum is largely converged at max- 
imum masses ~ 6 — 8 x 1O 14 /i _1 M0, with the exact values 
depending on the multipole - higher mass clusters contribute 
most to the lowest multipoles. The figure also shows the ratio 
of the standard deviation to the mean in the multipole bins 
/ = 2000 - 3000 and 3000 - 4000 as a function of the maxi- 
mum mass; for both multipole bins, the ratio increases as the 
maximum mass is increased. The skew also increases with 
the maximum mass of the cluster, approximately doubling 
between M max = 10 14 /i _:L M Q and 10 ls h^M®. Thus, it is 
the largest mass clusters that provide the largest amounts 
of non-Gaussianity within the power spectrum. 

3.6 Parameter dependence 

The five free parameters used in the cluster model (see equa- 
tions 12 and 15 1 are not yet constrained to high accuracy, 



with different simulations of galaxy clusters finding differ- 
ent values. Observations to measure the SZ effect are still in 
their early stages, and no large scale SZ surveys of clusters 
yielding many clusters have been carried out yet. To see the 
effect of the uncertainty in these parameters on our results, 
we vary the parameters over a wide range of possible values. 



3. 6. 1 Parameters of the Y-M relation 

The overall amplitude of the SZ effect is determined by Y t , 
which is fiducially set to 2 x 10 -6 Mpc 2 . As discussed earlier, 
Kay et al. (in prep.) find that this can be 1.9 or 2.3 x 10 -6 
Mpc 2 for non-radiative or preheated gas. Changing this pa- 
rameter simply scales the mean and standard deviation of 
the spectrum from the SZ effect, with the mean and stan- 
dard deviation being multiplied by (Y» :ncw / Fiducial) 2 - This 
is because the increase does not depend on the structure or 
distribution of the clusters. There is no effect on normalized 
quantities such as SBi/Bi and s from the SZ effect only. 
Due to the increase in the power from the SZ effect, though, 
the cross-over point between the CMB and SZ effect will 
be shifted to lower multipoles as Y* is increased, which will 
increase SBi/Bi and s at those lower multipoles. 

The evolution of the SZ effect as a function of mass 
is governed by 7. If the gas is assumed to only be gravita- 
tionally heated, then 7 = 5/3. This slope could be steeper, 
reflecting a reduced amount of gas mass in smaller clusters, 
or it could be flatter due to an increase in temperature in 
smaller clusters from extra, non-gravitational energy. Via 
X-ray observations of the inner region of galaxy clusters, 



Arnaud, Pointecouteau & Pratt (20071 find 7 = 1.82 



for the Yx-bA relationship in agreement with the preheat- 
ing result from Kay et al. (in prep.). As this is measured for 
the inner region, whilst the SZ effect extends much further 
out in the cluster atmosphere, this is likely to be an upper 
limit. As such, we use 7 = 1.5 and 1.8 as the expected range 
for this parameter. The effect of these values on the pow er 
spectrum and SBi/Bi is shown in the top row of Fig. 18 



4 Yx = M gas Tx is an X-ray proxy for the integrated SZ effect Y 
that has recently been shown to correlate well with M with little 



As 7 is increased, the peak in the SZ power spectrum 
shifts to lower multipoles, and SBi/Bi increases. This is be- 
cause the power from clusters with a mass over lO 14 /i _1 M0 
will be increased as 7 increases, whilst that from smaller 
mass clusters will decrease. The larger clusters are dominant 
at lower multipoles due to their greater angular size, so the 
peak of the SZ effect will naturally shift to lower multipoles. 
The effect is similar to simply increasing the mass of all the 
clusters. Section [5] showed that SBi/Bi increases as larger 
clusters are included in the maps; the same effect is present 
here, with high values of 7 providing the largest increase in 
the effective mass, hence the largest values of SBi/Bi. 

The SZ effect is expected to be self-similar, that is, it 
will depend on the redshift of the cluster via the evolution 
of the mean matter density within the virial radius of the 
cluster as well as the Friedmann equation. The parameter 
a models evolution beyond that from self-similarity, and is 
expected to be close to zero. Using numerical simulations, 



da Silva et al. (2004) find no evidence of evolution from self- 



similarity. 

To illustrate the effect of changing this parameter, we 
choose a modest range of values: a = —0.3 and 0.3; the 
results are presented in the second row of Fig. |18| Higher 
values of a will increase the SZ effect from clusters at the 
highest redshift. As the clusters at higher redshifts have a 
smaller angular size, this affects the power in the higher 
multipoles more than that in the lower ones. There are more 
clusters at higher redshifts, so when these provide increased 
power compared with the low redshift clusters the maps will 
become more similar to each other - it will become more 
difficult for a single cluster to stand out - hence the standard 
deviation of the maps will decrease. 

3.6.2 Parameters of the gas profile 

The final two parameters govern the profile of the SZ effect. 
The concentration of the gas within the cluster is determined 
by Cgas, and the slope of the cluster profile is governed by 



/3 (see equation 121. For (3 we use a fiducial value of 2/3, 



which is in agreement with observations of X-ray emission 
from nearby clusters by |Vikhlinin et al. (2006) that find 
(5 ~ 0.6—0.9 at rsoo (approximately half of the virial radius). 
Additionally, |Croston et al.| ( |2008| ) find a range of f3 ~ 0.37- 
0.81 for a range of radii between 0.3 and 0.8 -Rsoo- However, 
steeper values than these are expected in the outer regions 
where the gas pressure profile is expected to trace the NFW 



profile; this is demonstrated by Hallman et al. (20071, who 



scatter I Kravtsov, Vikhlinin fc Nagai|2006 1 



use simulations of clusters to compare X-ray and SZ profiles, 
finding that (3 = 0.88 for X-ray but 1.13 for the SZ effect. 
As a result, we investigate the effects of (3 = 1/2 and 1. 
For c gas , we use a fiducial value of 10 and also try 5 and 
20. Duffy et al. (in prep.) indicates that c gas is similar to 
com and depends on the cluster mass by around 10 per cent 
above 1O 14 M0, such that the values will lie around 5 — 10. 

The effects of c ga s and j3 are shown in the bottom two 
rows of Fig. |18| As these parameters only change the profile 
of the SZ effect, the average pixel value within the maps re- 
mains constant. The SZ effect is more concentrated within 
the centre of the cluster for higher values of /3 and c gas , with 
the effect of increasing the mean power spectra across all 
multipoles, but particularly at the higher multipoles, in a 
similar way to the effect of point sources on the power spec- 
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Figure 18. The effect on the power spectrum from changing the models of the parameters. From top to bottom the parameters are 7 
(the mass dependence) , a (the redshift dependence) , c gas (the gas concentration) and /3 (the outer gas profile) . The left hand plots show 
the mean power spectrum as a function of multipole I. The fiducial value is shown by the red solid line, with the blue dotted line showing 
the lower value and the black double-dotted line showing the higher value (see text for the parameter values). The data point is the range 
of values from the three values of ag in the 2000-3000 multipole bin. The right hand plots show SBi/Bi vs. the parameter values. The 
solid red line is for the multipole bin I = 2000 — 3000, and the blue dotted line is for I = 3000 — 4000. The dashed vertical lines show the 
three values used in the left-hand plots. The central vertical line is the fiducial value; the SBi/Bi for this value are between 0.23 — 0.27 
and 0.17 - 0.19 for the I = 2000 - 3000 and 3000 - 4000 bins respectively for the three values of erg. 
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trum (see section Eb. As the parameters are increased, the 
ratio of the standard deviation to the mean also increases, 
as the power from the central part of the largest (and rarest) 
clusters becomes more important at these multipoles. 

The assumption that all clusters are the same is un- 
likely to be true; in reality, there will be a certain amount 
of scatter within the cluster parameters, which may increase 
the standard deviation further when looking at smaller areas 
of the sky. Clusters are also likely to deviate from spherical 
symmetry. The extent of this scatter and its impact on the 
SZ effect is not yet known. 



3.6.3 Comparison to the effects of ag 

The changes in the gas physics have a comparable effect on 
the mean power spectrum as the range of o$ considered in 
this paper, especially for the outer gas profile /3 and the 
evolution of the Y — M relationship with redshift (a). At 
the same time, the gas physics can substantially increase the 
standard deviation of the power spectrum to a much larger 
degree than is possible with as, providing a much greater 
scatter in the power from different maps. These effects are 
obviously different to those from ag when the whole of the 
power spectrum can be considered, but will be difficult to 
separate when only a few bins are measured. 



4 POINT SOURCES 

Foreground emission from extragalactic point sources will 
contaminate the map. Synchrotron is the dominant process 
at frequencies below around 90 GHz, with emission from 
dust becoming important at higher frequencies. Although 
these sources are physically extended, this is generally much 
smaller than the beam size of the telescopes used to measure 
the CMB and SZ effect, so that they are well approximated 
as point sources. 

In our earlier results, we have not inserted point sources 
as CMB experiments generally attempt to subtract them. 
This procedure is, however, subject to uncertainties. Here, 
we discuss the effects of the inclusion of subdominant point 
source contributions for specific flux limits. We insert point 
sources statistically into these realizations, based on extrap- 
olation of their observed number counts. 

We calculate the total number of point sources with a 
flux density between 5 m i n and 5 ma x within an area AS! on 
the sky by 



iVtot = Af2 



dN 



dS v 



(21) 



where dN/dSv is the differential source counts as a function 
of flux density. To calculate the flux densities of individual 
sources, we randomly sample from a probability distribution 
defined by a normalized dN/dS v between the flux density 
limits. 

The flux densities can be converted to the brightness 
temperatures at a given frequency v using 



T V = S, 



./(- 



pixel 



dB\ 
df) 



(22) 



where #p ixol is the area of the pixel containing the point 
source and 
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Figure 19. Comparison of the 30GHz differential source counts 
from |Toffolatti et al.| < |1998[ ) rescaled by a factor of 0.7 (black 
dotted line) with the la ranges from measurements by (right-to- 
left) WMA P Ka band (Bennett et~ 
et al.||2002| blue), VSA (Cleary et al 



2003 



2005 



pink), DASI (Kovac 



et al.[2003| red) and GBT | |Mason et al.|2009] black) 



green), CBI \ Mason 



dB 
dT 



= 2k ( kTcms Y 

~ c 2 \ h J 



(23) 



is the differential of the Planck function with respect to tem- 
perature, where x is the dimensionless frequency, given by 
x — hpu/ksTcMB- 



4.1 Low frequency point sources 

For the low frequency point sources we use the 30 GHz model 
for dN/dS v from |Toffolatti et al.| ( p98l henceforth T98). 
This is an extrapolation of the observed number counts at 
1.4 GHz, which have been measured to much lower flux den- 
sities than the number counts at the frequencies of interest 
here. We have normalized the model by a factor of 0.7 (fol- 
lowing |clear^etr^n[2005| to bring it into closer agreement 
with measurements of the number count at 30 GHz (see Fig. 
191, although this still results in an over-prediction of the 
number of sources around 1 mjy compared with recent ob- 
servations using the Green Bank Telescope (Mason et al.| 
20091. 



To scale the flux densities to the desired frequency 
v, we use S^/Ssoghz = (^/SOGHz)". a is determined for 
each source using a Gaussian distribution about a mean 
a = —0.3, with o a = 0.36, for each source, which is the 
distribution for 15 to 30 GHz measured for the 9C sample 
(Cleary, private communication). We include point sources 
with fluxes between 10 -4 ' 5 and 10 5 mjy at 30 GHz; sources 
above this flux range are very rare. Sources with weaker 
fluxes are not well characterized and should have a negligi- 
ble contribution to the power spectra. 

In addition to the total number of point sources and 
their flux densities, we need positions. The simplest method 
is to distribute them randomly over the map, so that the 
probability of finding a source at a particular position, 
P(x,y), is constant across the map. However, galaxies emit- 
ting synchrotron radiation are thought to reside preferen- 
tially within galaxy clusters. Observations at 30 GHz by 



Coble et al. ( 2007 1 suggest that the number counts of these 
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Table 2. The statistics of the erg = 0.825 realizations containing CMB, SZ and point sources. The mean B; and the standard deviation 
SBi within the bin i are in £iK 2 ; the skew s is dimensionless. The top set are for 3x3 degree maps at 30 GHz using 0.2 mjy point 
sources; the middle for lxl degree maps also at 30 GHz and 0.2 mjy point sources, and the bottom set for 3x3 degree maps at 150 GHz 
with 2 mjy point sources. The three sets of values given are for no point sources (left), randomly positioned point sources (middle) 
and clustered point sources (right). Adding point sources increases the mean and standard deviation at high multipoles, as expected. 
Clustered point sources fill in some of the decrement from the SZ effect, reducing the mean and standard deviation. 



sources increase in the direction of known massive clus- 
ters; |Lin fc Mohr| ( |2007[ > come to a similar conclusion at 
1.4 GHz. To account for this, we correlate the positions of 
these sources with the galaxy clusters by using a power of 
the surface mass density T,(x, y) of the clusters. In particular 
we use the probability distribution 



P(x,y) = 



Y,{x,y) 



1/3 



(24) 



/ E(a;',y') 1/3 dx' dy 1 ' 

The specific choice of E(o;,j/) 1/ ' 3 brings the number counts 
towards clusters into reasonable agreement with those esti- 



mated by Coble et al. (2007) 



4.2 High frequency point sources 

Although negligible at the fiducial frequency of 30 GHz, 
dusty galaxies become important at higher frequencies. For 
these, we use the fit to SCUBA observations provided by 
Borys et al. (20031. This is in the form of a double power 



law, 

dN 
dSu 



So 



(ir+( 



(25) 



where the parameters at 350GHz are So = 1.8mJy, a = 1.0, 
p = 3.3 and N = 1.5 x 10 4 deg~ 2 . This is extrapolated to 



other frequencies using S v = v" 1 , where 7 = 2.5. We use the 
flux range 10 1 to 10 mjy at 350 GHz. These sources are 
expected to be clustered, but they are at higher redshifts 
than the galaxy clusters we are interested in. Thus, they are 
randomly distributed in these realizations. 



4.3 Power spectrum statistics with point sources 

Table [2] shows the power spectrum statistics for maps con- 
taining point sources below 0.2 mjy in addition to the CMB 
and SZ effect. The addition of point sources naturally in- 
creases the mean and standard deviation of the spectra at 
higher multipoles. The point sources also decrease the nor- 
malised skewness and the off-diagonal terms of the correla- 



tion matrix at 30 GHz (see Fig. 20 1 as they can reduce the 
amount of power from the SZ effect by filling in the decre- 
ments ( |Holder|"2 002). Clustering the point sources increases 
the chance of this happening, resulting in a slightly lower 
mean values for the power spectrum and a slightly reduced 
standard deviation. The effect at 150 GHz is similar to that 
at 30 GHz, except that the effect of clustering the low fre- 
quency point sources is reduced as the unclustered infrared 
sources are becoming more important. 

We do not attempt to account for the effects of imper- 
fectly removing strong point sources from the realizations, 
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Figure 20. The correlation matrix for 1000 realizations of the 
combined CMB, SZ effect and 0.2 mjy point sources at 30 GHz, 
using 3x3 degree maps. The addition of point sources decreases 
the off-diagonal terms (compare with Fig. |13| . 



which may result in a distribution of effective point sources 
with a flux below the cut-off flux. This could be an important 
issue, especially if the point sources lie in galaxy clusters, in 
which case the residual point source flux could fill in the 
decrement from the SZ effect. 



5 IMPLICATIONS FOR THE HIGH 
MULTIPOLE EXCESS 

The original, and highest significance, measurement of a pos- 
sible excess at large multipoles comes from the Cosmic Back- 
ground Imager. An excess was first announced by Mason et| 
al. ( 2003[| and was refin ed by the addition of more observa- 
tions by (Pearson et ah] ( |2003| ) and |Readhead et al.| ( |2004 



The final results were announced by Sievers et al.| (120091 
using the full 5 years of observations with the CBI. In to- 
tal, the CBI has observed five 5x5 degree fields and one 
5x6 degree field using mosaicked, shallow observations. Ad- 
ditionally, one 5 x 0.75 degree and three 0.75 x 0.75 degree 
"deep" fields were observed, which have reduced noise as 
they have been observed for longer. Following from |Bond et| 
al. (20051, Sievers et al. (2009) find that if this excess is due 



to the SZ effect, then erg must be between 0.9 and 1.0. 

We do not attempt to carry out a direct comparison 
between our results and the CBI measurements as we omit 
a number of factors that will be important in this compari- 
son. These include the difference in map sizes and the com- 
bination method - we investigate the statistics of 3 x 3 and 
0.75 x 0.75 degree maps separately, whereas the CBI mea- 
surement results from the combination of a number of differ- 
ent areas of the sky and different map sizes. We do not take 
into account the effects of the UV coverage, beam shapes, 
nor the window function for the observations. We also ex- 
clude any effects of noise, point source subtraction and any 
other foregrounds such as galactic emission or atmospheric 
effects. Finally, we only briefly cover the effects of gas physics 
and limits on the maximum cluster masses on the statistics. 
As such, our results are only indicative. 

We concentrate on two of the bins measured by the CBI, 
namely I = 2050 - 2350 and I = 2350 - 3900; these are the 
bins where the SZ effect will be most significant. For these 
bins, the CBI has measured 261±132^K 2 and 387±117/xK 2 




1 000 2000 3000 4000 5000 6000 7000 8000 90001 0000 
I 

Figure 21. The power spectra from the CMB+SZ effect for the 
three cosmologies, with the error bars showing the minimum and 
maximum values from all 1000 realizations with a map size of 3 X 3 
degrees. The black double-dotted points (offset by 100 multipoles 
for clarity) are for ag = 0.9, the red solid points for erg = 0.825 
and the blue dotted points (also offset by 100 multipoles) are for 
erg = 0.75. The two black data points are the latest values for the 
highest multipole bins from CBI jSievers et al.|2009| . The three 
solid lines are (from top to bottom) the mean power spectra from 
the SZ effect in with erg = 0.9 (black double-dotted line), 0.825 
(red solid line) and 0.7 (blue dotted line). 



respectively. They find similar values, but with larger error 
bars, when the deep fields are included or considered on their 
own, and find low significance excesses in each of the four 
fields seperately. The CBI data points are shown in Fig. |21| 
along with the power spectra from the SZ effect on its own 
for each three values of as, as well as the mean values from 
the combined CMB and SZ realizations for 3x3 degree 
maps with the spectrum binned with bin widths of 1000; 
the minimum and maximum values from the realizations 
are shown by the error bars. The highest CBI data point 
is clearly in excess over what would be expected from the 
mean and scatter of these realizations. 

To investigate the excess more systematically, we sam- 
ple the power spectra for the same bins as used by the 
CBI. The histograms of the power spectrum from the CMB 
and SZ effect within the multipole bin I = 2050 — 2350 for 
as = 0.825 and map sizes of 0.75 x 0.75 degree (representa- 
tive of the deep fields) and 3x3 degree are shown in Fig. 
|22| The measured value for the CBI power spectrum is com- 
patible with the simulations for this multipole range; 26.5 
per cent of the realizations have the same or higher value 
as the central value from the CBI measurement. The same 
histogram for the I = 2350 — 3900 bin is shown in Fig. |23| 
In this histogram, the distribution is significantly skewed 
towards higher values, however the probability of getting 
the CBI excess from the realizations is low, with only 0.36 
per cent of the realizations having the same or more power 
than the central CBI value. This would be reduced further 
if high-mass clusters are removed from the maps. 

Table [3] gives the values for the various components and 
maps sizes for the two CBI high-multipole bins. Although 
the mosaicked fields extend to 5 degrees across, we use 3x3 
fields due to the limits of the simulation box size at the 
highest redshifts. For the lower multipole bin, we find that 
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Table 3. Values for the mean and standard deviation at 30 GHz for the I = 2050 - 2350 (top) and I = 2350 - 3900 (bottom) CBI bins 
for 3x3 degree (1000 realizations) and 0.75 X 0.75 degree (16000 realizations) for various cosmologies and components, as well as the 
minimum and maximum values found and the normalized skew. The mean Bi and the standard deviation SBi within the bin i are in 
fiK 2 ; the skew s is dimcnsionlcss and and the values of S max given are in mjy. The components are the Cosmic Microwave Background 
(CMB), the SZ effect (SZ), point sources (PS) and clustered point sources (CPS). The CBI data points are 261 ± 132 and 387± 117 fj,K 2 
for the lower and higher bin respectively; the final column gives the percentage of realizations at or above the central values of these 
measurements. 
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Figure 22. The statistics of the combined CMB and SZ effect 
realizations with erg = 0.825 for the multipole bin 2050-2350. 
The large grey histogram is for the 0.75 X 0.75 degree realizations 
(bin width of 75 £iK 2 ); the small, darker histogram is for the 
3x3 degree realizations (bin width of 9 /iK 2 ). The CBI value 
for this bin is 261 ± 132/iK 2 (mean and ±lcr represented by the 
three dashed vertical lines); 26.5 per cent of the 0.75 X 0.75 degree 
realizations have this central power or higher, and 5 per cent have 
376.5^K 2 or higher. 

the central value can be easily obtained using the CMB and 
SZ effect with as = 0.9, and realizations remain with the 
same amount of power for as = 0.825. Due to the large 



Figure 23. As Fig. [22] but for the multipole bin 2350-3900 with 
bin widths of 28 and 3 /jK 2 for the lxl and 3x3 degree real- 
izations respectively. The CBI value for this bin is 387 ± 117 /jK 2 
(mean and ±lcr represented by the three dashed vertical lines); 
0.36 per cent of the 0.75 X 0.75 degree maps have this central 
power or higher, with 5 per cent having 137.7/xK 2 or higher. 

error on the measurement, however, it is in agreement with 
a large range of values of erg. For the higher multipole bin, 
there are no realizations with the central power for any of 
our three values of as using 3x3 degree realizations. There 
are 0.75 x 0.75 degree realizations that match the central 
value, even down to as = 0.75, but these are very unlikely, 
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making up less than a percentage of the realizations for the 
lower values of as . 

Adding point sources to the realizations at the level of 
0.2 or 1 mjy significantly increases the mean of the real- 
izations in both bins, greatly increasing the number of re- 
alizations with the central power or higher as measured by 
the CBI in both bins. Clustering the point sources, however, 
slightly decreases this probability. 

Due to the uncertainty in the parameters controlling 



the power and profile of the SZ effect (see Section 3.6 1, the 



mean power spectra given here are uncertain to around 50 
per cent. If this effect is in the positive direction, then this 
would significantly increase the number of realizations which 
agree with the CBI excess measurements. We change two of 
the parameter values, setting Y» = 2.3 x 10 -6 Mpc 2 and 
(3 — 1 in two different sets of realizations. We find that the 
change in Y* has little effect, however the change in f3 has 
a dramatic effect on the standard deviation and the skew of 
the distributions, as well as on the mean. 

As shown earlier in Fig. |I7| a large part of the skew and 
increased standard deviation compared to the mean comes 
from the largest clusters in the maps. Although there is no 
formal mass limit for the CBI fields, they were selected to 
avoid known large clusters, making it unlikely that any local, 
high mass clusters are present in the fields. This will have the 
effect of reducing the number of realizations matching the 
power measured by the CBI in accordance with the results 
of section . 

In addition to the CBI, the Berkeley-Illinois-Maryland 
Association (BIMA) interferometer operating at 28.5 GHz 



has also measured an excess. Dawson et al. (20061 surveyed 
eighteen 6.6 arcminute fields with a total area of ~ 0.2 
square degrees, and found 2201i2q/iK 2 with an average mul- 
tipole of I = 5237 and a FWHM of 2870. We model this as a 
constant bin between I = 3800 and 6670. For a thousand re- 
alizations with as = 0.825 and the same map size as BIMA's 
6.6 arcminute field of view, we find a mean within this bin 
of 41 fiK 2 , a standard deviation of 82 /iK 2 and a normalized 
skew of 5.5. Using just the mean and standard deviation, 
this makes the central value of the measurement a ~ 2.2 a 
excess. We find that 3.3 per cent of the realizations have 
greater than 220 /iK 2 . 

The Sunyaev-Zel'dovich Array (SZA) has also measured 
the power spectrum at higher multipoles, putting an upper 
limit of 149/^K 2 (95% confidence) on the power spectrum 
from secondary anisotropies between I = 2000 — 6000 ( Sharp 
et al. [ 2009) using fourty-three 12 arcminute fields. As 68 per 
cent of their multipole coverage is between I = 2929 — 5883, 
we use a constant bin with a width of I = 2900 — 5900. 
We find a mean of 50/iK 2 , with a standard deviation of 64 
and a skew of 5.8. 5.6 per cent of the realizations have a 
power greater than 149^tK 2 . This is compatible with the 
measurement from SZA. 

At higher frequencies, the Arcminute Cosmology 
Bolometer Array Receiver (ACBAR) has also measured an 
excess of 34±20^K 2 at 150GHz in the range I = 1950-3000; 
this 1.7 a excess is compatible with the CBI and BIMA 
measurements if the signal has the spectral distribution of 



For comparison, however, we find that the SZ effect alone 
contributes 9.6/^K 2 with a standard deviation of 2.2^iK 2 and 
a skew of 0.94 within the bin I = 1950 - 3000 for 3 x 3 de- 
gree maps. We note that QUaD ( Friedman et al.|2009 l also 
measure the 150 GHz power spectrum at high multipoles; 
they find values that are compatible with ACDM on its own 
with no SZ effect required. 



6 CONCLUSIONS 

We have simulated the microwave sky, including the CMB, 
SZ effect and point sources, and have used PlNOCCHlO to 
generate large numbers of cluster catalogues with realistic 
distributions for three different values of erg. Using these 
maps, we have investigated the statistics of the power spec- 
trum between multipoles of 1000 and 10000. 

We find that the inclusion of the SZ effect increases 
the standard deviation of the power spectrum by a factor 
of 3 over that expected from cosmic variance, in agreement 
with the predictions from an analytical calculation based on 
the halo formalism. The mean and standard deviation vary 
as l//gky as expected, and scale approximately as <r| over 
the range of values sampled here. We also find that the dis- 
tributions are non-Gaussian, and are skewed by large mass 
clusters, with the degree of this skewness increasing as the 
map size is decreased. Additionally, we find that correlations 
between galaxy clusters play a small role in the statistics of 
the power spectrum at the level of ~10 per cent. 

Several instruments have measured an excess at high 
multipoles, which may be due to the SZ effect with a large 
value of as- We cannot explain the central values of these 
measurements with the range of as investigated here, how- 
ever the increased standard deviation and the presence of 
skewness in the distribution means that these measurements 
could be explained by a lower value of as than has been 
suggested so far. There is also a large uncertainty in the pa- 
rameters describing the cluster gas physics, which can have 
a large effect on the mean of the distributions, comparable 
to that from the different values of as, and can also signif- 
icantly effect the standard deviation and the skew of the 
distributions. 

The next generation of CMB instruments are currently 
being commissioned, and are expected to provide more data 
at multipoles comparable to those probed by the CBI. The 
Planck spacecraft will measure the power spectrum from 
the whole sky out to multipoles of 2500 within the next few 
years, and instruments such as the South Pole Telescope 
Ruhl et al.|2004[), the Arcminute Microkelvin Imager 



Zwart et al.|2008|) a nd the Atacama Cosmology Tele- 



the SZ effect (Reichardt et al. 20081. These measurements 



cover nearly 700 square degrees over 10 fields; we are limited 
by the map size of the Pinocchio simulations, which are 9 
square degrees, so cannot make a prediction for these fields. 



(SPT; 
(AMI 

scope (ACT; Fowler et al.||2004 1 will observe large numbers 
of galaxy clusters using the SZ effect. These measurements 
will provide much more information on the SZ effect and 
may provide a resolution to the discrepancy in as from the 
measurements to date. Our results should also be relevant 
to these observations. 
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APPENDIX A: ANALYTICAL CALCULATION 
OF THE SZ POWER SPECTRUM 

Rather than calculating the mean power spectrum from 
maps of the SZ effect, it can be evaluated directly from the 
number density and cluster profile using the halo formalism 
(see for example Cooray| 2001 Komatsu & Seljak 20021 un- 
der the assumption that the clusters are Poisson distributed. 
The mean binned spectra can be calculated by 



igbin 



1(1 + 1) 

2tv 



dV dn o i „ , i / a , n 
dz dM yt ' v ' 
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where i is the bin number and the sum is over the multipoles 
I in the bin. 

dv( z) c 3 (/; E( z ' )dz f 

dz E(z)H$ 1 ' 

is the comoving volume of the Universe at a given redshift, 
where E(z) is as in equation [4] and dn/dM is the comoving 
cluster number density. The parameter yi is the value of the 
Fourier transform of the cluster model at the multipole I, 
given by 

yi(M,z) = 2Tvg(x) / <p c y(<p c )Jo((l + 0.5)tp c )dip c , (A3) 
Jo 

where g(x) determines the frequency dependence of the SZ 



effect, as per equation 13 y(<p c ) is from equation 14 and Jo 
is zeroth-order cylindrical Bessel function. 

Using the covariance matrix, the standard deviation can 
also be calculated in a similar way. Here, the angular trispec- 
trum Tui is used, defined by ( |Cooray|2001| |Komatsu fc Sel-| 
jak|2002 ) 



"■rc^"** <A4) 

This is then combined with the expected Gaussian cosmic 
variance from the mean spectrum to calculate the covariance 
matrix Mu* via 

n r 1 / 2C;C// Tui \ f . 

Mu -7^{W^)Ai Su ' + ^)' (A5) 

where Al is the width of the multipole bin and / s k y is the 
fraction of the sky being considered. The standard devia- 
tion can then be calculated from the diagonal of the co- 
variance matrix, that is, ai = 1(1 + Mu/2-k . This effec- 
tively includes both Gaussian and Poissonian terms. Note 



that Zhang & Sheth (20071 add an additional component 
to the calculation of the variance to model the clustering of 
galaxy clusters but find that the Poissonian term dominates 
over the clustering term; for this reason we do not consider 
it here. 
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